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We present deep optical imaging of the z = 4.1 radio galaxy TN J 1338-1 942 
obtained using the Advanced Camera for Surveys (ACS) on-board the Hubble Space 
Telescope as well as ground-based near-infrared imaging data from ESO/VLT. The 
radio galaxy is known to reside within a large galaxy overdensity (both in physical 
extent and density contrast). There is good evidence that this 'protocluster' region is 
the progenitor of a present-day rich galaxy cluster. TN J1338 is the dominant galaxy in 
the protocluster, in terms of size and luminosity (in both the optical and near-infrared) 
and therefore seems destined to evolve into the brightest cluster galaxy. 

The high spatial-resolution ACS images reveal several kpc-scale features within 
and around the radio galaxy. The continuum light is aligned with the radio axis and 
is resolved into two clumps in the 2775 and 2:350 bands. These components have lumi- 
nosities ~ lO 9 I/ and sizes of a few kpc. The estimated nebular continuum, scattered 
light, synchrotron and inverse-Compton scattering contributions to the aligned contin- 
uum light are only a few percent of the observed total, indicating that the observed 
flux is likely dominated by forming stars. The estimated star-formation rate for the 
whole radio galaxy is ~ 2OOM yr _1 . A simple model in which the jet has triggered 
star-formation in these continuum knots is consistent with the available data. 

A striking, but small, linear feature is evident in the z 85 o aligned light and may be 
indicative of a large-scale shock associated with the advance of the radio jet. The rest 
of the aligned light also seems morphologically consistent with star-formation induced 
by shocks associated with the radio source, as seen in other high-2 radio galaxies (e.g., 
4C 41.17). An unusual feature is seen in Lyman-a emission. A wedge-shaped exten- 
sion emanates from the radio galaxy perpendicularly to the radio axis. This 'wedge' 
naturally connects to the surrounding, asymmetric, large-scale (~ 100 kpc) Lyman-a 
halo. We posit that the wedge is a starburst-driven superwind, associated with the first 
major epoch of formation of the brightest cluster galaxy. The shock and wedge are 
examples of feedback processes due to both AGN and star-formation in the earliest 
stages of massive galaxy formation. 

Subject headings: galaxies: active — galaxies: halos — galaxies: high-redshift — 
galaxies: individual (TNJ1338-1942) 



1. Introduction 

The most massive galaxies in the local Universe reside in the centers of rich clusters. Within 
the context of hierarchical models of biased galaxy formation, the mass of a galaxy and its clus- 
tering properties are naturally connected via the initial density fluctuations (e.g., Kaiser 1984). 
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Therefore, not only locally, but throughout cosmic time, massive galaxies mark the densest regions 
of the Universe. The study of young overdensities at high redshift ('protoclusters') then also traces 
the history of the future brightest cluster galaxies. 

Many observing programs, spanning wavelengths from radio to X-ray, have been devoted to 
identifying galaxy overdensities over a large range of redshift (e.g., Postman et al. 1996; Scharf 
et al. 1997; Stanford et al. 1997; Rosati et al. 1998; Oke, Postman, & Lubin 1998; Rosati et al. 
1999; Holden et al. 1999; Kurk et al. 2000; Pentericci et al. 2000a; Holden et al. 2000; Donahue 
et al. 2001; Francis et al. 2001; Stanford et al. 2002; Postman et al. 2002; Donahue et al. 2002; 
Mullis et al. 2003; Miley et al. 2004). To date the most distant protoclusters have been found at 
z ~ 5 (Shimasaku et al. 2003; Venemans et al. 2004). Do these very young overdensities already 
contain a dominant, massive galaxy? 

There are several observational clues to the mass of a high-redshift galaxy. One is the ob- 
served X-band magnitude, which probes the rest-frame optical out to z ~ 4. Another is the 
presence of a high luminosity active nucleus (a supermassive accreting black-hole), which implies 
the existence of a large spheroidal host galaxy at least locally (Magorrian et al. 1998; Gebhardt 
et al. 2000; Ferrarese & Merritt 2000). High-redshift radio galaxies (HzRGs) are bright at K and 
harbor powerful nuclei (Jarvis et al. 2001a; De Breuck et al. 2002; Willott et al. 2003; Jarvis & 
McLure 2002). Therefore the fields surrounding HzRGs are important targets for studying the ear- 
liest examples of massive galaxies and clusters. Using a narrow-band Lyman-a imaging program 
Miley and collaborators discovered an overdensity of star-forming galaxies around all four radio 
galaxies observed to sufficient depth, out to z — 5.2 (Venemans et al. 2004). The resulting set of 
protoclusters is the subject of several on-going studies. This discovery also implies that the radio 
galaxies are the seeds of brightest cluster galaxies. 

Brightest cluster galaxies (BCGs) are the most massive galaxies known in the local Universe; 
having stellar masses in excess of 1O 12 M (J0rgensen, Franx, & Kjaergaard 1996; Bernardi et al. 
2003). The luminosities and sizes of brightest cluster galaxies are not drawn from the same distri- 
butions as the majority of the galaxy population (Oegerle & Hoessel 1991). BCGs at low redshift 
lie on the extrapolated Fundamental Plane of elliptical galaxies (Oegerle & Hoessel 1991). The 
surface-brightness profiles of some BCGs, the cD galaxies, extend out to hundreds of kiloparsecs. 
These shallow power-law stellar envelopes blur the distinction between the galaxy and the diffuse 
intracluster light. Such extreme sources are clearly very useful laboratories for studying the pro- 
cesses inherent to massive galaxy and cluster formation. In fact, several authors have shown that 
the observed build-up of BCGs can provide key constraints on the hierarchical theory of galaxy 
formation (Aragon-Salamanca, Baugh, & Kauffmann 1998; Burke, Collins, & Mann 2000). The 
present discrepancies between the predicted (using semi-analytic models) and observed abundance 
of massive galaxies imply that fundamental processes are not being accounted for in the current 
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models (Cole et al. 2000; Baugh et al. 2003; Somerville et al. 2004). One possibility to solve 
this discrepancy is to postulate the existence of strong interactions between accreting black-holes, 
star-formation and their host galaxy and surroundings (i.e., "feedback"). 

Galaxies which host powerful radio sources are peculiar in several respects. The most striking 
property is the radio-optical 'alignment effect'. This effect is the strong tendency of the rest- frame 
ultraviolet continuum light of the radio host to be aligned with the axis defined by the radio source. 
Several explanations for this behavior have been put forward: recent star-formation induced by 
the radio jet, nebular continuum from emission line gas which is photo or shock ionized by the 
AGN, light from the central engine scattered into the line of sight by either dust or electrons, and 
inverse-Compton scattering of the microwave background or other local photon fields (McCarthy 
et al. 1987; Chambers, Miley, & van Breugel 1987; Daly 1992; Dickson et al. 1995). These expla- 
nations of the alignment, but particularly the jet-induced star- formation, are excellent examples of 
feedback on the galaxy formation process. 

The giant (~100 kpc) Lyman-a emitting halos surrounding distant radio galaxies may be an 
observable consequence of feedback from galaxy formation (e.g., van Ojik et al. 1997; Reuland 
et al. 2003). The enormous line luminosities of these objects, often in excess of 10 44 ergs s -1 , 
imply they are massive reservoirs of gas (Kurk et al. 2000; Steidel et al. 2000). What is the 
origin of this gas; outflow from the galaxy or infall of primordial material? There are several 
pieces of circumstantial evidence that these halos are connected with the AGN. The halos are 
often aligned with the FRII radio axis, and in some cases Lya emission is directly associated with 
radio structures (Kurk et al. 2003). Perhaps photoionization by the AGN or shocks due to the 
radio source expansion are responsible for the extended line emission. Spectroscopically, Lya 
absorption is seen in addition to the bright halo emission (e.g., van Ojik et al. 1997; Wilman et al. 
2004). There is a good correlation between the amount of absorption and the size of the radio 
source (Jarvis et al. 2001b). A possible scenario is that a neutral hydrogen shell initially surrounds 
the radio source, but is subsequently ionized during the growth of the radio source. There are also 
some halos known which do not contain a bright radio source (Steidel et al. 2000) suggesting that 
perhaps the halo phenomenon is associated with the more general processes of galaxy formation, 
rather than being specific to active nuclei. 

The radio galaxy TN J1338-1942 (z = 4.1; De Breuck et al. 1999) resides in one of the 
youngest protoclusters known (Venemans et al. 2002; Miley et al. 2004). This galaxy lies within 
a large Lya halo which shows unusually asymmetric morphology when compared to other similar 
radio sources (Venemans et al. 2002). In this paper we present high spatial-resolution HST Ad- 
vanced Camera for Surveys (ACS) imaging of this radio galaxy. These data provide the clearest 
view of a young brightest cluster galaxy to date. Images in four broadband filters ((7475 ,r 62 5 ,2775 
and zs5o ) have been obtained. The resulting magnitudes and colors have been used to apply the 
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'Lyman-break' technique to select galaxies at the same redshift as the radio galaxy (Miley et al. 
2004). The exquisite spatial resolution of HSTI ACS allows us to study the detailed morphology of 
the radio galaxy. Using these data we present a scenario which describes the observed morphol- 
ogy (both continuum and Lya) and the measured kiloparsec-scale colors and magnitudes within a 
self-consistent formation framework for TN J1338— 1942 (hereafter TN J1338). 

This paper is structured as follows: we describe the observations and data reduction in § 2, 
we present the results of our analysis of the combined multi-wavelength dataset in § 3, and discuss 
these results in § 4. We adopt the 'concordance' cosmology (Spergel et al. 2003) with Vt m = 0.27, 
= 0.73 and H = 71 km s _1 Mpc -1 . Within this cosmology the angular scale at the redshift 
of TN J1338, z = 4.1, is 7.0 kpc arcsec -1 . We use the AB magnitude system (Oke & Gunn 1983) 
except where noted. 



2. Observations and Data Reduction 
2.1. ACS Imaging 

The ACS data of TN J1338 were taken with two primary goals in mind: first, to enable 
color- selection of faint protocluster members using the 'Lyman-break' technique and second, to 
investigate the detailed morphological structure of the brightest protocluster galaxies, including the 
radio galaxy. To achieve these goals, images were taken in 4 broad-band filters: #475 (F475W), r 62 5 
(F625W), i 775 (F775W) and z 850 (F805LP). The <7 475 -band is below the Lyman-break for galaxies 
at z > 4.1. The r e2 5 -band contains (and for the case of the radio galaxy, is dominated by) the 
Lyman-a emission line. Both the i 775 and z 850 bands are relatively unaffected by bright emission 
lines (see Table 1). However, we note that one emission line, CIII]A1909, may be affecting the 
2 850 -band morphology of the radio galaxy. The #475 , r e2 5 and i 775 observations were carried out 
between 8 and 12 July 2002 and the z 850 images between 1 1 and 12 July 2003 with the Wide Field 
Channel of the ACS. The total observing time of 18 orbits was split between the broad-band filters 
9400 seconds in g 475 , 9400 s in r 625 , 1 1700 s in i 775 , and 1 1800 sec in z 850 . Each orbit was split 
into two 1200 sec exposures to facilitate the removal of cosmic rays. g 475 dropouts were selected 
using the g 475 , r 62 5 and i 775 bands as described in Miley et al. (2004) and Overzier et al., in prep. 
This effort was very successful and confirmed the presence of a galaxy overdensity around TN 
J1338. In this paper we present the first discussion of the radio galaxy itself and its role as the 
dominant protocluster galaxy. 

The ACS data were reduced using the ACS pipeline science investigation software (Apsis; 
Blakeslee et al. 2003a), developed for the ACS Guaranteed Time Observation (GTO) program. 
After the initial flat-fielding of the images through CALACS at STScI, the Apsis processing steps 
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include the empirical determination of image offsets and rotations, the rejection of cosmic rays, 
the combining of images through drizzling, and object detection and photometry using SExtractor 
(Bertin & Arnouts 1996). The final images have a scale of 0''05 per pixel, and (2a) limiting (AB) 
magnitudes of 28.46 (g 475 ), 28.23 (r 62 5 ), 28.07 (i 775 ), and 27.73 (z 850 ) in 0.2-arcsec 2 apertures 
(corrected for Galactic extinction). 



2.2. VLT Optical Imaging and Spectroscopy 

Deep VLT/FORS2 images (6f8 x 6'.8) of the TN J1338 field were made in March 2001 in 
both the broad-band R filter (A c = 6550 A, AA F whm = 1650 A) and a custom narrow-band filter 
(A c = 6195 A, AAfwhm = 60 A) to target the redshifted Lyman-a emission line (Venemans 
et al. 2002). The la limiting surface-brightnesses are 28.6 and 29.2 per arcsec 2 for the narrow 
and broad bands respectively. These images were used to identify candidate Lya emitting galaxies 
(Venemans et al. 2002). For the current paper these images are used to elucidate the larger scale 
structure of the Lyman-a emission and to compare it to the structures seen in the ACS images. 

A follow-up spectroscopy program using FORS2 in multiobject mode was carried out in May 
2001. These spectra have a dispersion of 1.32 A pixel -1 (using the 600RI grism) and cover the 
wavelength range from 5300 to 8000 A. Candidate Lya galaxies were placed on two slit-masks 
(Venemans et al. 2002). The radio galaxy itself was included on both slitmasks providing a deep 
spectrum along the radio axis which covers both the Lya and CIVA1549 emission lines. 



2.3. VLT NIR Imaging 

.fCj-band images of TN J1338 were obtained in two separate observing runs. One on March 
24-26 2002 collected 2.1 hours of total exposure time using ISAAC on UT1 of the Very Large 
Telescope (VLT). The second run, using the same instrument, was done in service mode at VLT 
between the nights of 27 May and 13 June 2004. The total exposure time for this run was 5.7 hours. 
All the data were processed, sky- subtracted and combined using the XDIMSUM package within 
IRAF 1 (Tody 1993). The final image has a scale of 0.148"per pixel, a seeing of 0"5, and a (2a, 
1 sq. arcsec aperture) limiting magnitude of 25.6. ISAAC has some geometrical distortion across 
the face of the detector. We have not corrected the distortion in detail because the K s morphology 
of the radio galaxy matches features seen in the continuum ACS observations. 



'IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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The fCs-band is the only band which probes wavelengths longward of the 4000 A break, 
beyond which an old stellar population should dominate the emergent flux (see Figure 1). It should 
be noted, however, that the bandpass is not entirely at A > 4000 A. However, these data still provide 
a crucial point on the spectral energy distribution (SED) of the radio galaxy. 



2.4. Radio Imaging 

The radio source TN J1338 was originally selected because of its ultra-steep spectrum (be- 
tween 365 MHz and 1 .4 GHz) which has been shown to be an indicator of high redshift (De Breuck 
et al. 2000b). The first radio data were culled from the Texas and NVSS catalogs (Douglas et al. 
1996; Condon et al. 1998). Follow-up observations were made with the VLA at 4.71 and 8.46 
GHz in March 1998 (Pentericci et al. 2000b). The noise levels are 25 and 50 /iJy beam -1 for the 8 
and 5 GHz maps, respectively. The resolution is 0!'23 for the 8 GHz and 0"43 for the 5 GHz map. 
These are the primary radio data used in this paper (see overlay in Figure 2). The radio source 
has three distinct components at both these frequencies; the northwest (-S^tqhz = 21.9 mJy) and 
southeast (£ff GHz = 1.1 mJy) lobes (separated by 5"5) and the likely radio core (S^q Hz = 0.3 
mJy) located very close (1"4) to the northwest lobe (De Breuck et al. 1999). The radio source is 
highly asymmetric, with the NW lobe nearly 20 times brighter at 4.7 GHz than the SE lobe. The 
radio asymmetry may indicate an asymmetry in the ambient medium (McCarthy, van Breugel, & 
Kapahi 1991). 



2.5. Image Registration 

The center position, orientation and angular-resolution all differ between the ground-based 
and space-based data sets. To facilitate comparison amongst these datasets we registered them all 
to a common reference frame and pixel scale. The ACS images were registered with respect to 
each other in the ACS GTO pipeline. We chose the ACS r 62 5 frame as this common grid. For 
the ground-based R and K s data this registration requires interpolation from plate-scales of 0"201 
and 0!'148 to 0!'05 per pixel respectively. The shifts, rotations and rebinning were all done in a 
single interpolation step using the IRAF tasks geomap and geotran. We fitted a general coordinate 
transformation using 2nd order polynomials in both axes. The rms deviations of the data from the 
fits were of order 0.3 input pixels (0"04) for the K s image and 0.07 input pixels (O('Ol) for the 
FORS2 narrow and broad band images. The calculated transformations were done using 'sine' 
interpolation within geotran. At least 15-20 unsaturated stars were matched within the ACS r 625 - 
band and each of the ground-based images to calculate the appropriate transformations. We applied 
these transformations to the ground-based R, narrow-band and K s images. 
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The optical/near-infrared frame is defined by stellar positions and that of the radio image is 
defined by the positions of radio point sources (quasars). The radio image of the TN J1338 field is 
sparse, so a direct matching of sources will not produce a robust transformation between the two 
frames; particularly when we do not want to use the radio galaxy itself. Therefore the accuracy 
limit to the radio-optical registration is determined solely by the systematic error in the optical 
reference frame. To better quantify this we have used two different optical frames, the USNO and 
the GSC-2.0, and compared the astrometric solutions using both to the radio data. This gives us 
a conservative amplitude of 0"3 to the error on the position of the radio core with respect to the 
optical structures in the ACS and ground-based data. Therefore we can confidently associate the 
radio core with the region near the peak of the K s -band flux and at one tip of the ACS galaxy. 



2.6. Continuum subtraction 



To model the continuum in the r 62 5 we have used a power-law extrapolation from the relatively 
emission-line free i 775 and z 850 bands. While these filters are not strictly line-free due to the pres- 
ence of CIVA1549 in the i 775 band, HeIIA1640 in both passbands and, to a lesser extent, CIII] A1909 
in the z 850 band, none of these lines are expected to dominate the continuum (based on spectro- 
scopic data, see Table 1). We assume that the continuum follows a simple power-law in F\(oc X 13 ) 
that extrapolates through to the ACS r 62 5 bandpass, VLT R and narrow-bands. We have accounted 
for the intergalactic absorption shortward of the emission line and the relative throughputs of the 
filter curves. Light at wavelengths shorter than Lyman-a is easily absorbed by neutral hydrogen 
located between the source and the observer. This will greatly affect the amount of continuum light 
detected in any bandpass shortward of the emission line. We have therefore adopted the model of 
the intergalactic hydrogen optical depth presented in Madau (1995); r cfr = 0.0036(^) 3 46 . This 
is the optical depth due to Lyman-a forest lines, and not due to higher order Lyman series or metal 
lines. These two opacity sources make negligible contributions to the total optical depth at this 
redshift (Madau 1995). We integrate this attenuation over the filter curve shortward of Lya to de- 
termine the amount of continuum flux absorbed in the intergalactic medium (IGM). This correction 
decreases the amount of continuum by 23.9%, 16.9% and 36.2% in the ACS r 62 5 , VLT R and VLT 
narrow-bands respectively. In addition, the continuum flux detected in the images is affected by 
the filter throughput curves. The extrapolation from continuum measured in some filters (in our 
case, the ACS i 775 and z 8 50 ) to other bandpasses depends on the relative filter curves. We have 
used the total integrated throughputs to correct for the differential sensitivities. The final value of 
P for the entire radio galaxy is —1.32 (or a = 0.62, where F v oc v a ). 
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3. Results 

The ACS and VLT images are shown in Figure 3. The rest- frame ultraviolet morphology of 
TN J1338 is complex and multi-faceted. The radio galaxy exhibits the usual alignment between 
its continuum, line emission and radio axis (Chambers, Miley, & van Breugel 1987; McCarthy 
et al. 1987; Best, Longair, & Rottgering 1998). Two kiloparsec- sized clumps along the radio axis 
dominate the continuum structure in the i 775 and z 850 bands. This is similar to many powerful 3CR 
radio galaxies at z ~ 1 (Best, Longair, & Rottgering 1998; Zirm 2003; Zirm, Dickinson, & Dey 
2003). In the r 62 5 band, the ACS image reveals both concentrated and diffuse Lyman-a emitting 
regions. In this section we quantify the line and continuum flux distributions for the radio galaxy 
by performing photometry in a set of varied apertures. 



3.1. Estimating the Contributions to the Aligned Light 

The presence of a powerful radio source has several effects on the emitted spectrum. As 
mentioned above, at least three of these effects tend to align the observed continuum emission with 
the axis defined by the double radio lobes. Any remaining continuum flux we attribute to a young 
stellar population. 

Ionized gas not only emits emission lines, but also continuum photons from two-photon re- 
combination, bremsstrahlung and standard recombination. For a given gas temperature and den- 
sity, this nebular continuum spectrum can be calculated (see Figure 1). Observationally, the nor- 
malization of the nebular spectrum is determined by using the emission line spectrum. For TN 
J1338— 1942 we have a spectrum which was taken along the radio axis and spatially-averages the 
entire optical extent of the radio galaxy (RG). Ideally one would measure Balmer recombination 
emission lines. These lines directly correspond to the nebular continuum. At z = 4.1, these lines 
have shifted out of the optical window visible from the ground. As a 'Balmer proxy' we use the 
HeIIA1640 emission line (?), which has a flux of « 1.7 x 1(T 16 ergs s" 1 cm~ 2 (2004, C. De 
Breuck, private communication). Assuming the HeII/H/5 ratio (= 3.18) from a high-redshift radio 
galaxy composite spectrum (McCarthy 1993), we estimate the H(3 flux to be 5.4 x 1CT 17 ergs s _1 
cm~ 2 . If these emission lines arise in a 15,000 K gas, at low enough densities where collisional 
de-excitation is negligible, then the corresponding nebular continuum flux densities in the r 62 5 , 
i 775 and z 850 bands are 1.5 x 1CT 31 , 3.1 x 1CT 31 and 3.8 x 10~ 31 ergs s _1 cm -2 Hz 1 respectively. 
So, at its brightest the nebular continuum averaged over the entire galaxy is only 2; 850 (AB) =27.4, 
much fainter than even the individual components of TN J1338. We subtract the nebular continuum 
from all quoted magnitudes and fluxes by scaling the subtracted amounts for the sub-components 
by their estimated Lya flux. 
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The same population of electrons that is responsible for the radio synchrotron emission can 
also inverse-Compton (IC) scatter ambient photon fields. The primary seed photons are those 
making up the cosmic microwave background (CMB), with secondary contributions from the syn- 
chrotron photons themselves (synchrotron self-Compton; SSC) and other AGN emission. By as- 
suming equipartition to calculate the magnetic field in the northern radio lobe of TN J1338, we 
find a P-field on the order of a few hundred micro-Gauss. Using the extrapolation to the rest-UV 
as calculated by Daly (1992), this translates to an IC-CMB contribution of only a few xl0~ 31 
ergs s _1 cm~ 2 Hz -1 , or fainter than 28 magnitudes (AB). The estimated energy density in the syn- 
chrotron field is similar to that of the CMB at this high-redshift, so SSC adds a similar amount to 
the aligned continuum. These two contributions combined are negligible. As with the spectroscopy 
above, these values are spatially-averaged over the entire galaxy. 

Optical spectropolarimetry of the entire galaxy reveals a polarization of 5 ± 3 % (2004, C. De 
Breuck, private communication). This value is much lower (about a factor of two) than similar 
measurements of z ~ 1 RGs (e.g., Dey et al. 1996; Cimatti et al. 1996, 1997; Solorzano-Inarrea 
et al. 2004), but similar to other powerful radio galaxies at z ~ 4 (Dey et al. 1997). The amount of 
scattered light is related to the percent polarization by the intrinsic polarization, p. p depends on 
the type of scatterer (either dust or electrons) and the geometry of the scattering, neither of which 
is well-constrained for TN J 1338. Therefore we can only put the lower limit on the percentage 
of scattered light of 5%. However, we also note that for z ~ 1 3CR radio galaxies the observed 
polarizations are high and suggest that P is also high. We therefore conclude that for TN J1338 
the aligned continuum contains some scattered light (~ 10%), but is substantially diluted by unpo- 
larized sources. We have already shown that neither the nebular continuum nor the IC scattering 
can account for this dilution. 

Clear evidence for the presence of young stellar populations in radio galaxies has proven dif- 
ficult to find. The same spectral region, the ultraviolet, where massive stars are brightest coincides 
with the bright region of the AGN spectral energy distribution. Therefore, the detection depends 
on being able to find stellar- specific features in very deep spectra (e.g., the SVA1502 photospheric 
absorption feature in 4C 41.17 at z = 3.8 Dey et al. 1997). Existing spectra of TN J1338 do not 
show such features. In this case the presence of young massive stars must be inferred from the UV 
excess after subtracting the other known Contibutors to the aligned light. As we have shown above, 
there seems to be such an excess in TN J 1338, which implies the existence of many young stars. 
We examine this result further below. 
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3.2. Large Aperture Photometry 

To perform integrated photometry over the whole galaxy we used a "Kron" aperture (Kron 
1980) to approximate a total galaxy magnitude in the ACS and VLT images. This enables easy 
comparisons between the space and ground-based observations of TN J1338. This aperture is op- 
timized to be as large as possible while still retaining a high signal-to-noise ratio. The magnitudes 
were determined using SExtractor (Bertin & Arnouts 1996). The Apsis pipeline "deblends" the 
radio galaxy into two distinct objects. We therefore reset the SExtractor deblending parameters 
to maintain the radio galaxy as a single object. We overplot the aperture in Figure 3. The 'total' 
magnitudes within this aperture are listed in Table 1 . 

The radio galaxy is > 1.4 magnitudes brighter than the next brightest "dropout" or Lyman-a 
emitting galaxy (within similarly defined apertures). This is true even in the K s band where the 
light from older stars is presumed to dominate the emission rather than processes related to the 
AGN. Therefore the radio galaxy is securely identified as the dominant, and likely most massive, 
galaxy within the protocluster and the probable progenitor of a present-day brightest cluster galaxy. 



3.3. Rest-frame Optical Surface-brightness Profile 

The K s -band samples the rest-frame optical continuum emission, mostly longward of the 
4000 A break, and shortward of the potentially bright [OIII]AA4959, 5007 and H(3 lines. At these 
wavelengths it is likely that the galaxy luminosity is dominated by older, low-mass stars. We show 
the final ISAAC image of the radio galaxy in panel (e) of Figure 3 and the raw surface-brightness 
profile of the radio galaxy and a star in this band in Figure 4. This profile is centered on the peak 
of the K s -band light. It is clear that the galaxy is resolved, but it is not possible to measure the size 
of the galaxy robustly from these data. A simple r^-law fit roughly estimates the effective radius 
at 0'.'7 or ~ 5 kpc. This is considerably smaller than the corresponding radii of local BCGs (e.g., 
Graham et al. 1996). In fact, the observed morphology in the K s image is very similar to that seen 
in the ACS i 775 -band. The galaxy is highly elongated and is nearly perfectly aligned with the radio 
axis. There are no bright emission lines within the bandpass, and none of the continuum processes 
discussed above are bright in K s . Therefore, we conclude that the stellar distribution is aligned 
with the radio axis. 
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3.4. Multi-component Decomposition and Diffuse Light 

The rest-frame ultraviolet morphology of TN J1338 can be decomposed into several discrete 
components. To better understand the morphology of the galaxy we have extracted multi-band 
photometry of these individual components. We divided the radio galaxy into distinct regions 
both in the line-dominated r 625 -band and in the continuum-dominated i 775 -band. The resulting 
segmentation map is shown in Figure 5. The regions are numbered 1 thru 6. Regions 1, 2 and 3 
('wedge') are found only in the emission-line dominated r 62 5 -band. Conversely, regions 4 and and 
5 are most prominent in the continuum bands (i 775 and z$ 50 ). Finally, region 6 contains a linear 
feature in z 850 . 

For each region we have measured the magnitude in each ACS band and in the ACS Lya 
(continuum subtracted r 62 5 ) image. Using these values we have derived the UV continuum slope 
and line flux independently for each portion of the galaxy. In addition, we have measured the 
magnitudes in the fT^-band for the two primary continuum clumps visible in the ACS data (regions 
4 and 5). The (i 775 — K s )ab color for region 5 is ~ 0.2 magnitudes redder than region 4. The sum 
of the light contained within these clumps is about one magnitude fainter in the i 775 -band than 
the total radio galaxy, indicating that there is some "diffuse" light detected in these high angular- 
resolution images. This remains the case in the zq^q -band with the extended light apparently 
having the same color as the mean color of the entire galaxy. The magnitudes and Lya fluxes are 
presented in Columns 3-6 of Table 2. 



3.5. Azimuthal Binning 

While the aligned light is a common feature of high-redshift radio galaxies, the wedge of ex- 
tended line-emission to the southwest of the radio galaxy (region 3 in Figure 5) is unusual. We have 
measured the azimuthal surface-brightness distribution of this feature by extracting photometry in 
angular bins. The bins are shown in Figure 6. The angular regions are sized and placed to cover 
the entire visible extent of the wedge and to be small enough to accentuate the internal structure of 
the feature. The flux in each bin was summed using the IRAF task 'polyphot' . The corresponding 
errors were calculated by performing the same aperture photometry on the error maps. The pixels 
within the radio galaxy (as defined by the extent of the i 775 -band continuum) were masked out. 
The resulting flux histogram is shown in Figure 7 for the wedge (solid line), the opposite side of 
the galaxy in the line image ('anti-wedge', dashed line) and for the i 775 -band continuum in the 
wedge region (dotted line). 

The bins covering the wedge show a clear excess of flux with respect to both the opposite 
side of the galaxy and to the continuum. There also appears to be some structure to the azimuthal 
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profile. The azimuthal profile has relatively sharp-cutoffs at either side. There is no evidence for 
limb-brightening toward the edges of the wedge. However, the spatial resolution (~ 1.5 kpc/bin) 
is insufficient to resolve a thin shell of limb-brightened emission. 



3.6. Semi-circular Annuli 

To measure the radial profile of the wedge we have constructed semi-circular annuli which 
cover the same region(s) as the angular bins above (Figure 8). We have increased the width of the 
annuli towards the outer edge of the wedge to maintain constant signal-to-noise ratio per annulus. 
The annular photometry again shows that the line emission is significantly brighter in the wedge 
than on either the opposite side of the galaxy or in the continuum. The surface-brightness profiles 
for the wedge (solid line), 'anti-wedge' (dashed line) and continuum (dotted line) are shown in 
Figure 9. The errors for each data point were derived by performing identical photometry on the 
corresponding error map. There does not seem to be evidence for limb-brightening toward the 
outer edge of the wedge. Also shown are power-law profiles of the form SB(r) oc r~ a with a 
values of 1 and 2 (dot-dash lines). The wedge profile has a profile much closer to the a = 1 curve. 
We shall consider these profiles further in the discussion. 

4. Discussion 

The unparalleled spatial resolution provided by HST and the Advanced Camera for Surveys 
has allowed us to observe kpc-scale structures within the radio galaxy TN J1338— 1942 at z = 
4.1. The rest-frame ultraviolet continuum and line emission of the host galaxy is morphologically 
complex and consists of several distinct components. In this section we discuss these features, 
their implications for the formation of this galaxy and attempt to construct a possible scenario for 
the on-going processes in this source. 

4.1. On-going Star-Formation 

Are we witnessing the formation of the bulk of the final stellar mass in this radio galaxy? 
Previous studies based on ground-based optical and millimeter imaging suggest that the current 
star-formation rate (SFR) within TN J1338 is very high, of order several hundred solar masses 
per year (Venemans et al. 2002; De Breuck et al. 2004). The current study not only confirms this 
estimate, under different assumptions, but also determines the spatial gradients in star-formation 
activity. This new spatial information allows us to constrain both the timescale and possible phys- 
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ical origin of the current star-formation episode. 

There are two methods to estimate the SFR in TN J 133 8 using the ACS imaging data; one 
using the Lya luminosity and the other using the UV continuum luminosity. Both estimates are 
based on the assumption that the UV light we see from the galaxy is due solely to star-formation. 
As we have shown earlier, the UV continuum is relatively unpolluted by the non-stellar contributors 
to the aligned light. Conversely, for the line emission, the high observed equivalent-widths argue 
in favor of non-stellar excitation (e.g., Chariot & Fall 1993). Therefore, the line emission provides 
an upper limit (modulo dust extinction) to the SFR, while the UV continuum provides a lower 
limit. The comparison of these two derived values provides an estimate of the mean SFR for each 
portion of the galaxy. 

For a given Lya flux we calculate the corresponding Balmer line flux and hence the SFR 
(Kennicutt 1998). This SFR is a lower limit locally due to Lya being a resonantly scattered line 
and very susceptible to dust. For example, the Lya flux of the wedge (region 3 in Figure 5) is 
3.3 x 1CT 16 ergs s _1 cm~ 2 corresponding to an intrinsic luminosity of 5.2 x 10 43 ergs s _1 at z — 4.1. 
If we adopt a Case B Lya/Ha ratio of 8.7 (Brocklehurst 1971), and use the relation between Ha 
line flux and SFR, we find that region 3 is forming stars at a rate of 40 M Q yr _1 . The analogous 
results for the other regions are listed in Table 2. The total "Kron" Lya SFR is 290 M yr _1 . The 
sum of the individual SFRs does not equal this total; underpredicting it by m 90 M yr -1 . There 
is a considerable amount of low surface-brightness line emission even within the "Kron" apeture 
(~ 18 x 7 kpc). We have not corrected for the observed HI absorption (De Breuck et al. 1999; 
Wilman et al. 2004). 

The 2 775 and z 850 band probe the rest-frame UV of TN J1338-1942 at ~ 1500 A and ~ 1775 
A, respectively. We can use the ACS continuum magnitudes to estimate the current star formation 
rate, under the assumptions that the UV luminosity is dominated by the light from late-O/early-B 
stars on the main sequence. We measure total ("Kron") star formation rates of 86 and 96 M Q yr _1 
for the i 775 and z 850 bands, assuming a Salpeter IMF. The derived SFRs are lower limits, since 
they are dependent on the amount and distribution of dust present in the UV emitting regions. The 
slope of the UV continuum can also be used to measure the extinction. We measure the slope of 
the continuum from the i 775 — z S50 color and use a template spectrum of a typical star-forming 
galaxy redshifted to z — 4.1 to convert the measured slope to a color excess, E(B — V). For 
the template spectrum we have used the stellar population synthesis models of Bruzual & Chariot 
(2003) to create a typical "Lyman-break" galaxy spectrum, having an exponentially-declining star 
formation history (with time constant, r = 10 Myr), an age of 70 Myr, 0.2 Z & metallicity, and a 
Salpeter initial mass function. The parameters of this template are taken from the best-fit SED at 
z ~ 3 of Papovich et al. (2001). We varied the dust content by applying the attenuation curve of 
Calzetti et al. (2000) to this template. We find E(B - V) = 0.12, yielding a dust-corrected SFR 
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of ~ 220 M & yr _1 in good agreement with the emission line estimate above. If we repeat this 
calculation for each discrete region of the galaxy (allowing for E(B — V) and SFR to change for 
each; see Table 2), we again find evidence for diffuse UV light and star-formation. The sum of the 
SFRs for all the regions again falls short of the total by a factor of 4.5, "missing" 170 M yr _1 . 



4.1.1. Shocks and Jet-Induced Star-Formation 

One of the primary explanations for the alignment effect is that the passage of the radio jet 
through the interstellar gas induces star-formation (e.g., Rees 1989). Strong large-scale shocks 
associated with the expanding radio source overpressure molecular gas clouds which then collapse 
to form stars. The presence of powerful shocks in radio galaxies at z < 2 has been inferred via 
their ultraviolet emission line ratios (e.g., Best, Rottgering, & Longair 2000; De Breuck et al. 
2000a). For TN J1338 most of the important diagnostic emission lines are unobservable from the 
ground. However, we can use morphological information to search for possible signatures of shock 
processes in this radio galaxy. 

In these respects, useful analogies can be drawn between TN J1338 at z = 4.1 and 4C 41.17 
at z — 3.8. The shock properties of 4C 41.17 were studied and modeled in detail by Bicknell et al. 
(2000). In addition to HST imaging of 4C 41.17, these authors also used deep emission line spectra 
and high angular-resolution radio imaging to study the relationship between the radio source, the 
gas and the stars. We apply a similar analysis to TN J1338. 

A bow-shock is formed at the tip of the advancing radio jet (cf. Figure 3 in Bicknell et al. 
(2000)). Due to the unresolved radio structure we cannot robustly determine the location of the jet 
interaction. We assume that region 4 in the continuum ACS image is the primary site where the jet 
has shocked or is still impacting the gas. This continuum knot is very blue and has a morphological 
structure which is suggestive of jet-cloud interaction; namely, a paraboloid oriented along the radio 
axis. From the spatially-resolved optical images we estimate the interaction area to be ~ 2 x 10 44 
cm 2 by assuming the emission we see is emitted by a spherical shell. The jet is most likely well- 
collimated and therefore the area of the jet itself is much smaller than the total; we assume 10%. 
If we assume that most of the momentum flux of the jet is dissipated in this interaction, the shock 
velocity (see Eqn. 1 of Bicknell et al.) will be greater than: 

v sh > 300 - 2000F^ 4 2 6 /3 j - 1/2 n- 1/2 km s" 1 (1) 

F E)46 is the energy flux of the jet in units of 10 46 ergs s _1 , /?j et is the relativistic Doppler parameter, 
and n H is the hydrogen density per cm 3 in the cloud. Comparison between model and observed 
CIVAA1548, 1550 doublet fluxes can help to constrain the pre-shock gas density and the energy 
flux of the jet. Assuming that the ACS UV/Lya image also the spatial distribution of CIV, we can 
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constrain the area from which the line is being emitted. We use the observed CIV/Lya flux ratio to 
convert the Lya image to a "CIV" image. If we follow Bicknell et al. (2000) and take A p to be the 
projection of the true area of the shock A sh in region 4 and predict the CIV line luminosity for our 
estimated shock velocity we find: 

L(CIV) « 2 x 10 42 [^p]^(^)- 1 ergs (2) 

a(CIV) is the radiative efficiency of the CIV doublet. By comparing this to the observed line 
luminosity (4.8 x 10 42 ergs s _1 ) we estimate that the pre-shock electron number density to be on 
the order of ?t.h = 3 cm~ 3 . 

We have estimated the star-formation rate produced by the jet-cloud interaction in region 4. 
To determine whether these shocks can lead to star-formation rates in excess of a couple dozen 
solar masses per year (for region 4) we parameterize the jet-induced star-formation as follows. 

SFR jct = 26(^ T )(^)(^)A sM4 ^ h , 1O ooM yr" 1 (3) 

We have taken the ambient pre-shock gas number density, p, to be 3 cm~ 3 as inferred above using 
the line luminosity. The area of the shock front is taken from the image itself and the shock velocity 
was assumed to be between the extremes possible in equation 1 above. The efficiency with which 
the shocked gas is converted into stars is denoted by e and the gas volume filling factor by f gas . 
The assumption of order unity filling factor is most likely incorrect globally in the galaxy, but more 
likely to be realistic in this limited region where the gas is approximately uniform. We conclude 
that jet-cloud interactions could be responsible for the observed star-formation in TN J1338. 

The existing radio imaging data of TN J 133 8 does not have sufficient angular-resolution to 
make direct correspondence with specific features in the ACS images. However, since the radio 
structure clearly overlaps the optical galaxy the comparison with 4C 41.17 and our assumption of 
a physical connection between the radio and optical structures are likely to be justified. 



4.2. The 'Wedge' 

The coupling and regulation of star-formation and nuclear activity in forming galaxies and 
their interaction with their environments (i.e., "feedback") is a key issue in modern astrophysics. 
At low redshift, imaging and spatially-resolved spectrocopic observations of starbursting galaxies 
show that supernovae-driven outflows, so-called "superwinds", are common among this population 
(e.g., Heckman, Armus, & Miley 1990). Whether such outflows actually achieve escape velocity 
and enrich the inter-galactic medium (IGM) is still an outstanding question (Heckman 2000; Heck- 
man et al. 2000; Martin 2004). However, these winds, along with nuclear outflows, are the only 
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processes observed to transport material into the outer halos of galaxies and are therefore prime 
candidates for injecting the metals and energy that are observed in the IGM. At high-redshift, 
where these processes are likely to be even more prevalent due to the higher global star-formation 
rates, the strongest evidence for the presence of outflows is spectroscopic. However, the spectro- 
scopic features observed are of ambiguous origin and could be due to inflow, outflow or rotation 
(e.g., van Ojik et al. (1997) but see also Adelberger et al. (2003)). Spatially-resolved imaging of 
the emission-line gas can provide a more certain indication of outflow, if the gas is collimated 
or exhibits the bi-polar morphology of low-redshift superwinds. In the next sections we consider 
several possible origins for the wedge: the photoionization cone of an AGN or young stellar pop- 
ulation, in-situ star-formation, scattering by dust, or an ionized outflow associated with a starburst 
(i.e., a superwind). 



4.2.1. Photoionization 

In several low-redshift Seyfert galaxies cone-shaped regions of high ionization are observed, 
consistent with photoexcitation by an active nucleus (e.g., Wilson & Tsvetanov 1994). In some 
cases these cones extend to distances of 15-20 kpc from the nucleus, similar to the size of the Lya 
wedge seen in our ACS image (Wilson & Tsvetanov 1994). Both the ionized cone and our wedge 
have high equivalent-width. We derive a lower-limit to the Lya rest-frame equivalent-width of 
650A for the wedge emission. Furthermore, powerful radio galaxies are known to photoionize their 
surroundings based on emission line diagnostics and imaging of low-redshift sources. Generally 
both the line and UV continuum emission are elongated and aligned with the radio axis, particularly 
at z > 0.7 (McCarthy et al. 1987; Chambers, Miley, & van Breugel 1987). 

Conversely however, the principal wedge axis is perpendicular to that of the radio source. The 
unified model for AGN posits that the observed radiation is anisotropic (Antonucci 1993) due to 
a combination of obscuration close to the nucleus and the intrinsically anisotropic radiation. For 
radio-loud galaxies this preferential radiation axis is traced by the line connecting the dual radio 
lobes. The misalignment of the wedge therefore argues against photoionization due to the AGN 
or shocks due to the radio jet. Furthermore, the most likely position of the accreting black hole 
powering the radio emission and therefore also the primary source of hard ionizing radiation is 
where the radio core and K s -band surface-brightness peak coincide. The apex of the wedge does 
not coincide with this position. It is possible that a second AGN (this one radio-quiet), coinciding 
with region 4 in Figure 5, could ionize the wedge. However, this additional black hole would have 
to have its primary axis roughly perpendicular to the radio-loud AGN and be much less luminous 
at K s -band. While not impossible, this explanation is ad hoc and not preferred. 
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4.2.2. In Situ Star Formation or Galaxy Interaction 

On-going star-formation within the wedge itself, and perhaps extending into the outer Lya 
halo (outside the "Kron" radius), would also produce bright line emission. However, there is a 
robust upper-limit to the Lyman-a equivalent width produced by normal massive stars of W\ = 400 
A (Chariot & Fall 1993). For the wedge we estimate a significantly higher equivalent width. We 
also note that dust extinction will decrease the observed equivalent width from its true value; the 
resonant scattering of Lya photons increases their optical depth relative to continuum photons. 
This equivalent width argument also applies to tidal debris ejected from the galaxy via a merger or 
interaction. It therefore seems unlikely that stars are directly responsible for the wedge emission. 



4.2.3. "Superwind" : Comparison with M82 

In Heckman, Armus, & Miley (1990) the authors use their observations of local starburst 
galaxies to determine the minimal condition for driving a galactic- scale outflow powered by super- 
novae explosions, or a 'superwind'. They phrase this criterium in terms of the star-formation rate 
per unit area (EsfrX and emprically determine the minimum to be Esfr > 0.1 M yr" 1 kpc~ 2 . If 
we adopt this minimal value for TN J 133 8 and apply it to the resolved area where we see the wedge 
emerging from the galaxy (assuming a circular region seen in projection), we derive a lower-limit 
to the SFR over this same area of 1.5 M Q yr -1 . Above we have shown that the SFR for this galaxy 
greatly exceeds this limit even only in region 4 where the wedge may originate. 

Galactic-scale winds have been observed in detail around local starburst galaxies, a well- 
studied example being M82 (e.g., Heckman et al. 1990). Morphologically, these superwinds are 
bipolar structures emanating from the galaxy nucleus and along the minor axis of the galaxy. They 
are detected as emission-line filaments, extended X-ray lobes and bipolar thermal dust emission 
(Heckman et al. 1990). In M82, the emission-line gas is photoionized in the innermost regions and 
primarily shock-excited in the outskirts. The optically emitting gas flows from M82 in filaments 
which trace the biconical surface. The outflows have double-peaked emission lines (due to the two 
surfaces of the cone being separated in velocity) and also, depending on geometry and spectral 
resolution, emission lines with blue-shifted absorption. On the basis of these spectroscopic sig- 
natures, galactic-winds seem to be a generic feature of high-redshift star-forming galaxies (Pettini 
et al. 2001). All the existing spectroscopy of TN J1338 has been taken along the radio axis where 
the dynamics of the gas are presumably dominated by the AGN outflow and jet-cloud interactions 
(see above and e.g., Villar-Martin et al. 1999; Solorzano-Inarrea, Tadhunter, & Axon 2001). So 
we must rely on morphology alone to infer the presence of a superwind emanating from this radio 
galaxy along the wedge. 
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Comparison of the recent ultraviolet GALEX image of M82 (Hoopes et al. 2005) with our TN 
J1338 image reveals a striking degree of similarity (Figure 10). The scales of the two outflows are 
somewhat different; in M82 the narrowest collimated section is only 1.5 kpc across, while in TN 
J1338 the similarly defined region is approximately twice that. The surface-brightness profile is 
also similar between the wedge and the M82 far-UV outflow (see Figure 9). The shallow drop-off 
of emission-line surface-brightness is consistent with having shock ionization dominate at these 
large radii rather than photoionization by a central source. 

Why do we see only one side of the presumably bipolar outflow? This could be due to obscu- 
ration of the line emission on one side of the galaxy. If the radio galaxy is flattened (perpendicular 
to the plane of the sky) and inclined with respect to the line-of-sight then any dust in the galaxy 
would naturally obscure the side tilted away from the observer. This is the primary explanation 
for the observed asymmetry in the bipolar superwind in M82 (Shopbell & Bland-Hawthorn 1998). 
Alternatively, the lack of line emission on the northeast side may be due to a strong gradient in the 
ambient gas density. 

Alternatively, if the pressure is higher on the northeast side of the galaxy, a situation which 
could arise due to the motion of the radio galaxy within the protocluster medium, i.e., ram pressure, 
then the outflow would be impeded on this side of the galaxy. Marcolini, Brighenti, & D'Ercole 
(2004) have made simulations (albeit for dwarf galaxies) that show that motion through the IGM 
does not greatly affect the dynamics of galactic outflows until the ram pressure becomes compara- 
ble to the static thermal pressure of the galactic ISM. This result should be extendable to the case 
of TN J1338. To estimate the thermal pressure in the galaxy we must first estimate a density for the 
gas. If we assume that the observed hya emission is due to Case B recombination at T = 15000K, 
then we can use the fiducial Ly a/Ha ratio of ~ 10 to deduce the number of ionized hydrogen 
atoms. Given the observed geometry of the wedge we can assume further that the emitting gas is 
contained in the surface of a cone with half-angle 30° and length 10 kpc. The thickness of the gas 
layer cannot be larger than a few hundred parsec due to the absence of significant limb-brightening 
in the azimuthal profile of the wedge. Consequently, the number density and mass of ionized 
hydrogen are: 

n e = l-OL^^y-^cm- 3 (4) 

M H „ = 7.6 X 10V^^,4lKonc, k pc3 M (5) 

Where /jl p is the mean particle mass, which we have taken to be the proton mass, L H/ 3,4i is the 
H/3 luminosity in units of 10 41 ergs s _1 . The true electron density is likely to be higher than this, 
but only within smaller clouds or filaments as is seen in the M82 outflow. The resulting thermal 
pressure is equivalent to the ram pressure produced by a relative velocity of 300 km s _1 (if the 
surrounding gas has a density l/1000th of that in the wedge). This is some evidence that the radio 
galaxy is not in the center of the galaxy overdensity, and that motion towards the center would 
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provide the requisite ram pressure (Interna et al. 2005, in prep.). We conclude that the wedge 
is likely to be a supernovae-driven ouflow, with the current episode of star-formation possibly 
triggered by the radio jet. This superwind is one-sided due to ram pressure inhibiting the flow on 
one side. 



4.3. The Outer Lya Halo 

In Figure 1 1, the extended Lya (out to ~ 100 kpc), as detected in the very deep VLT narrow- 
band image, is shown as contours. The TN J1338 halo is somewhat has an asymmetric 'plume' 
which is aligned with the radio axis. It is clear from the underlying ACS image showing the wedge 
that there is a natural connection between the high surface-brightness wedge (out to 20 -30 kpc) 
and the larger scale lower surface-brightness halo along the southwest direction. However, the halo 
at larger distances appears to be aligned with the radio axis of TN J1338. What is the relation, if 
any, between the wedge and the large-scale Lya structure? 

If the wedge is an outflow, as we conclude above, the resulting bubble will stall at some 
radius where gravity and the amount of swept up intergalactic matter balance the input energy. 
The gas deposited at this radius will naturally flow along the boundary of the excavated cavity and 
follow any density gradients in the ambient medium. The halo-radio alignment would then be a 
natural consequence if lower density regions were preferentially along the radio axis. This would 
be the case if either the current radio source extended further out in radius than our current radio 
observations indicate, or that the radio source was previously (either during this same accretion 
episode or during an earlier one) much larger and had excavated the region along the radio axis. 
However, there is no evidence in the current radio data for a relic radio source at larger distances. 
In any case, if the ionized gas has originated in the starburst, the observed alignment implies that 
the AGN had already 'imprinted' the region before the starburst was triggered. We discuss this 
possibility and its implications further in the next section. 

4.4. A Self- Consistent Scenario 

The host galaxy of the powerful radio source TN J1338— 1942 at z = 4.1 is unique. It is 
arguably the youngest brightest cluster galaxy known to date, and has been the subject of several 
multi-wavelength investigations. In this paper we have presented imaging from HSTI ACS which 
reveals several interesting morphological and broad-band spectral features in this radio galaxy. In 
this section we attempt to construct a plausible and self-consistent story of the past, present and 
future of TN J 133 8. 
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The host galaxy of TN J1338 appears to be forming stars at a high rate. None of the non-stellar 
processes known to produce the alignment effect in other galaxies can be dominant in this case. The 
ACS data presented here reveals a morphology that is consistent with most of the star-formation 
being triggered by the passage of the radio jet. In Figure 12 we show a color-color diagram for the 
discrete regions within the radio galaxy. The color difference between regions 4 and 5 may be an 
age effect (region 4 is also bluer in i 775 — K s ). We use the overplotted model colors in Fig. 12 
to derive an age difference between the two regions. The model has constant star-formation with 
0.4 solar metallicity and E(B — V) — 0.1, the ages are labeled at timesteps of 1, 10, 100 and 
1000 Myr. We estimate from the comparison between the model points and the data that the age 
difference is between 25 and 200 Myr. This matches within the (large) errors the shock travel time 
from the radio core (in region 5) to the jet-cloud interaction in region 4 (a distance of ~ 7 kpc 
with a 300 km/s shock). We do not see very extended star-formation which may be triggered by 
the expanding radio cocoon. The energy injection from the AGN may be rather isolated to these 
few nodes along the radio jet itself. Therefore, we suspect that the large-scale Lyman-a gas has an 
origin apart from the AGN, namely in the newly-formed stars. 

Once the prodigious star-formation is initiated, the supernovae explosions expel the ionized 
gas into the surrounding media by means of a superwind. In addition, ram pressure is stripping 
the gas from the star-forming regions as the radio galaxy moves through the ICM/IGM. This gas 
remains ionized primarily via shock excitation at the interface between the outflow and the ambient 
intergalactic gas. At a distance of approximately 20-30 kpc from the nucleus, the outflow reaches 
pressure equilibrium with the IGM. The gas rapidly cools at this interface, consistent with the 
boundary to the Lya halo observed at this region. This cooling gas will follow the density gradients 
in the IGM. The very low surface-brightness line emission seen to the northeast of the radio galaxy, 
which is aligned with the radio axis, is probably gas which originated in the outflow, cooled, and 
is now being re-excited by low-luminosity shocks associated with the radio source. 

The observed rest-frame B magnitude of TN J1338 is approximately 1 .5-2 magnitudes brighter 
than that of the six brightest cluster galaxies at z ~ 1 observed with ACS (?). If we use the same 
star-formation model as above (constant for 1 Gyr and 0.4 solar metallicity) and age the galaxy 
from 1 Gyr at z = 4.1 to 5.2 Gyr at z = 1 the galaxy fades by 3.3 magnitudes. This would imply 
that some additional star formation or merging must occur during those 4.2 Gyrs for TN J1338 to 
match the luminosity of the z ~ 1 BCGs. This is certainly not surprising. It is interesting to note 
as well, that several of these BCGs have a nearby bright companion, of almost equal luminosity 
(AM < 0.1 mag), with which it seems to be destined to merge. The resulting increase in lumi- 
nosity would nearly make up the 1 magnitude of "extra" fading seen for the 1 Gyr constant star 
formation model. 

It is debatable whether one should attempt to draw conclusions for an entire population of 
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sources (either BCGs or radio galaxies in this case) based on observations of a single example. 
TN J1338 may be a galaxy which in a special phase of its evolution, or alternatively, it may be a 
special source whose history cannot be generalized to describe other galaxies. However, studies of 
ensembles of radio galaxies, including their luminosity functions and duty-cycles, suggest that the 
space density of radio source hosts at high-,2 are roughly in agreement with the density of BCGs at 
low redshift and the density of non-RG overdensities at z ~ 3 (?Venemans et al. 2002) 

5. Conclusions and Future 

The host galaxy of powerful radio source TN J 1338-1 942 shows signatures of several feed- 
back processes which connect the black-hole, stellar host and intergalactic medium. The elongated 
and multi-component ultraviolet continuum is aligned with the FRII radio axis and is likely to be 
due to emission from young stars being formed along the jet axis. Interpretation of this light in 
terms of jet-induced star- formation is consistent with the observations. There is, however, also 
evidence for considerable star- formation outside the highest surface-brightness regions. If the cur- 
rent star-formation rate has been constant over the jet travel-time from the radio core to the site 
of the presumed jet-cloud interaction, this process could have produced in excess of 10 11 M of 
stars. Data from the Spitzer Space Telescope will allow us to determine the total stellar mass of TN 
J 1338 and help verify our hypothetical star-formation history. 

We interpret the ACS wedge of Lya emission as a superwind driven by the winds and super- 
novae explosions associated with prodigious star-formation activity. This outflow connects with 
the larger-scale Lya halo. The initial source of the halo gas is then within the starburst, and is 
possibly enriched. Deep spectroscopy from the ground along the wedge axis would help us de- 
termine the ionization mechanism more definitively. In particular, covering the Lya, CIV, Hell 
and CIII] emission-lines may enable us to also measure the enrichment of the outflowing gas. An 
on-going ACS program using the narrow-band (ramp) filter (PI: W. van Breugel) to image several 
high-redshift radio galaxies in Lya will discover how prevalent such wedge features are in this 
population. 
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Table 1 . "Kron" Aperture Photometry 



Bandpass 


AB Magnitude 


Estimated Line 






Contamination 21 


3475 


25.92 ±0.28 


< 0.05 mag a 


»*625 


22.46 ±0.01 


~ 1.3 mag (Ly«) 


*775 


23.23 ±0.03 


~ 0.3 mag (CIVA1549,HeIIA1640) 


^850 


23.11 ±0.04 


~ 0.2 mag (HenA1640,Cm]A1909) 


K S 


21.9 ±0.2 


< 0.01 mag b 



a Calculated using the observed spectrum of TN J1338-1942 
b Calculated assuming the composite spectrum of McCarthy (1993) 
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Table 2. Six Easy Pieces 



Number 


Piece 


AB Magnitude 


Lyman-a Flux 


Lyman-a Rest-frame 


Extinction 


Star Formation Rate 






»-625 a 


*775 


Z850 


(10~ 1S ergs/s/cm 2 ) 


Equivalent Width (A) 


E(B - V) 


UV Continuum (Lya) b 


1 


Line 1 


27.14 


26.82 


26.80 


1.17 


448 


0.02 


3(17) 


2 


Line 2 


27.66 


27.34 


27.32 


0.80 


492 


0.02 


2(12) 


3 


Wedge 


26.19 


26.36 


26.81 


4.20 


645 


0.0 


3(61) 


4 


Continuum 2 


25.74 


25.48 


25.51 


3.73 


390 


0.0 


9(54) 


5 


Continuum 1 


25.33 


24.98 


24.92 


1.69 


121 


0.06 


24 (24) 


6 


Shock 


26.54 


26.20 


26.15 


2.45 


536 


0.04 


7(35) 



a Calculated for the Lya subtracted image, where the line flux is determined from power-law continuum fits to the 4775 and 2:350 bands for each 
piece 

b The star formation rate as calculated by the UV continuum magnitude and the Lya flux (in parenthesis) 
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Fig. 1. — Rest-frame wavelengths covered by the ACS and ground-based near-infrared filters. 
Shown left to right are the filter curves for the F475W, F625W, F775W, F850LP and K s band- 
passes. Also shown for comparison are the characteristic spectra of the nebular continuum emis- 
sion (green), high-redshift radio galaxies in general (blue; McCarthy 1993) and an older (1 Gyr) 
stellar population (red). 
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Fig. 2. — Color composite image of the r 62 5 , £775 and z 850 images from ACS, with the VLA 5 GHz 
radio map overlaid. Notice the blue 'wedge' emanating from the south-western side of the radio 
galaxy, this appears only in the r 62 5 band and is likely due entirely to Lyman-a emission. The rest 
of the emission is clearly aligned with the axis defined by the two radio lobes. The continuum in 
the 2775 and z 850 bands is dominated by the two clumps along the same axis. 
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Fig. 3. — Four optical and one near-infrared images used in this paper. North is up and East to 
the left. Each cutout is 4"5 on a side. From left to right: the r 62 5 ACS image, the same image 
but continuum- subtracted, the 2775 ACS image, the z 850 ACS image and finally the Ks image from 
VLT. All the ACS images have been smoothed with a 1 .5 pixel Gaussian kernel. The arrow in panel 
(d) marks the linear feature which is likely a large-scale shock, which is also shown unsmoothed 
in the inset. The two lines indicate the radio axis and the ellipse is the Kron aperture used for 
photometry of the entire radio galaxy. 
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Fig. 4. — fTs— band surface-brightness profile for TN J1338 extracted in circular apertures (open 
circles) along with that of a star measured in the same manner (stars). The radio galaxy is clearly 
extended and there is some indication that the profile may follow a de Vaucouleurs law. However, 
these data does not have sufficient spatial resolution to decisively determine the form of the surface- 
brightness profile. 
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Fig. 6. — Azimuthal bins used to measure the profile shown in Figure 7 overplotted on the 
continuum-subtracted r 6 25 (left) and i 775 (right) images. 
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Fig. 7. — Azimuthal profile of the wedge (solid line), 'anti-wedge' (dashed line) and ^775 contin- 
uum (dotted line). Note the sharp cutoff on either side of the wedge and the possible sub-structure 
(radial 'filaments') in the profile. 
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Fig. 8. — Semi-circular annuli used to measure the radial profile of the wedge in Lya (left). The 
same annuli were also used on the 'anti-wedge' side of the galaxy, and on the i-band continuum 
image (right). North is up and East to the left. Each cutout is 5" on a side, with major tickmarks 
separated by 1 arcsec. The radial dependence of all three are shown in Figure 9. 
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Fig. 9. — Radial surface-brightness profile for the wedge (solid), 'anti-wedge' (dashed) and con- 
tinuum (dotted) annuli. The profile derived from the GALEX data for M82 normalized to the 1338 
wedge is plotted as a dot-dash line. The two straight dot-dash lines are power-laws of slope -1 and 
-2. There is a clear excess of flux in the wedge between radii of 3.5 and 10 kpc. The slope follows 
the GALEX profile and the a = — 1 power-law well out to 10 kpc. At 10 kpc, the wedge profile 
falls off rather sharply, but the errors are large. 
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Fig. 10. — left: Color composite of r 62 5 , 2775 and z 850 ACS images, right: GALEX near (red) 
and far-UV (blue) composite image of M82, showing the bipolar outflow. In this case the far-UV 
is interpreted as being continuum light from the starburst that is scattered by dust mixed with the 
outflowing gas (Hoopes et al. 2005). Note the great similarity between the morphologies of the 
two objects which suggests that the wedge emanating from TN J1338 is also a starburst-driven 
outflow. 
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Fig. 11. — Continuum- subtracted ACS r-band image overlaid with the VLT narrow-band image 
(contours). North is up and East to the left. The fiducial bar in the lower-right of the image is 
10 kpc long. The field-of-view is 15" x 22"5, each major tickmark is one arcsecond. Note the 
correspondence between the wedge and the larger scale structure of the halo seen in the ground- 
based image. The 'filament' extending to the north is aligned with the radio axis, although it 
extends far beyond the northern lobe. 
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Fig. 12. — r 62 5 - 2775 vs. i 775 - ^850 color-color diagram for the individual regions in TN J1338. We 
have excluded the wedge (region 3) from this plot since it is a clear outlier, r 625 - i 775 = —0.11 and 
i 775 - z 850 = —0.39. The estimated Lyman-a flux has been subtracted from the T625 magnitude for 
each region using a continuum fit to the i 775 - z 8 5o color. The overplotted red line shows the colors 
for a constant star-formation model with ages indicated in Gyr, E(B — V) =0.1 and Z = 0AZ Q . 



